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Abstract 

This paper illustrates the application of group theory to a quantum-mechanical 
three-dimensional quartic anharmonic oscillator with Oh symmetry. It is 
shown that group theory predicts the degeneracy of the energy levels and 
facilitates the application of perturbation theory and the Rayleigh-Ritz vari¬ 
ational method as well as the interpretation of the results in terms of the 
symmetry of the solutions . We show how to obtain suitable symmetry- 
adapted basis sets. 
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1. Introduction 


Quantum-mechanical anharmonic oscillators have proved useful for the 
analysis of the vibration-rot at ion spectra of polyatomic molecules 1-3]. Sev¬ 
eral aspects of such spectra as well as other molecular properties have been 
modelled by means of simple coupled oscillators mainly with cubic and quar- 
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tic anharmonicities[f - 201. These and other applications of the ubiquitous 
quantum-mechanical anharmonic oscillators motivated their study and the 
development of suitable methods for the treatment of the corresponding 
eigenvalue equations 21^70]. The list of papers just mentioned is far from 


being exhaustive and its purpose is merely to give an idea of the interest 
aroused by the anharmonic oscillators along the years. 

Some authors have taken into account the symmetry of the anharmonic 


oscillators in order to simplify their treatment 


59 


60 


and others re¬ 


sorted to the more formal point group symmetry (PGS) to classify the vibra¬ 


tional states 


15 


36 


37 


49 


59 


61]. Despite the great amount of information 


provided by PGS none of those papers exhibits a full application of such 
mathematical tool, except the analysis of two-dimensional cubic and quartic 

42], These pap ers motivated a re¬ 


oscillators by Pullen and Edmonds 


41 


cent application of PGS to a variety of Hermitian 


71 


anharmonic oscillators with space-time symmetry 


75 


741 and non-Hermitian 


79|. In the latter case 


PGS proved suitable for determining the conditions that complex anharmonic 
potentials should satisfy in order to support real eigenvalues. 

The aim of this paper is to reinforce the idea that PGS is most important 
in the study of quantum-mechanical anharmonic oscillators because we deem 


that such technique has been surprisingly overloo 


for the few works already mentioned above 


41 


ced or disregarded (except 


42 


71 


ia|)- 


I 11 Section [2] we introduce the model, a quartic anharmonic oscillator 
with considerably large symmetry described by the point group Oh- This 
problem was treated before by means of perturbation theory and a less 
formal approach to symmetry based on parity and coordinate-permutation 
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operations 59]. Here we classify the eigenstates of the unperturbed Hamilto¬ 
nian according to the irreducible representations (irreps) of that point group 
and predict the rupture of the degeneracy by the quartic perturbation. In 
Section [3] we apply perturbation theory through second order to verify the 
splitting of the eigenspaces predicted in the preceding section. In Section [I] 
we discuss the application of the Rayleigh-Ritz variational method with basis 
sets adapted to the symmetry of the problem. In three consecutive subsec¬ 
tions we discuss the harmonic oscillator basis set, the Krylov space and a 
non-orthogonal basis set. We show results illustrating the splitting of the 
unperturbed eigenspaces due to the quartic perturbation that breaks the 
symmetry of the system. In Section [5] we summarize the main results of 
the paper and draw conclusions. Finally, in the |Appendix A| we outline 


the main features of group theory that are necessary for the analysis of the 
present problem. 

2. Model 

Among the many models mentioned above we have chosen a three-dimensional 


quartic anharmonic oscillator already studied earlier by Turbiner 


59] 


H = p 2 x + p 2 y + p 2 z + x 2 + y 2 + z 2 + A [/3 (x 4 + y 4 + z 4 ) + x 2 y 2 + x 2 z 2 + y 2 z 2 ] , 
A > 0, (3 > 0, (1) 


who recognized that it exhibits the symmetry of a cube. The author stated 
that his classification based on parity and permutation operators was in¬ 
complete. In this paper we describe the symmetry properties of this model 


by means of the point group Oh 


8CH82]. In other words, this Hamiltonian 
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operator is invariant with respect to the symmetry operations indicated in 


the table of characters shown in Table |T] described in the Appendix A A 
detailed discussion of the construction of the matrix representation of the 
symmetry operations for the Oh point group is available elsewhere 73]. 


The eigenvalues E^ n n and eigenfunctions <fkmn{ x , U-, z ) °f H 0 = H (A = 0) 


are 


E 


kmn = 2(k + m + n)+3 


Vkmn{ x ,y,z) = <l>k{x)<l> m {y)<l> n {z), k,m,n = 0,1,..., 


( 2 ) 


where 4>j(q) is an eigenfunction of the one-dimensional harmonic oscillator 
Hho = P 2 q + Q 2 ■ Every energy level is T+ 1 )A+-) degenerate, where u = 
k + m + n. 

Throughout this paper we resort to the following notation for the permu¬ 
tation of a set of three real numbers 


{a,a,a}p = {a, a, a} 

{a,b,b} P = {{a,b,b},{b,a,b},{b,b,a}} 

{a, b , c} P = {{a, b, cj, {c, a, 6}, {b, c, a}, {b, a, c}, {c, b, a}, {a, c, b}}, (3) 
where it is assumed that o ^ 6, a ^ c, and b ^ c. It enables us to express 
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the symmetry of the unperturbed eigenfunctions (A = 0) as 


{2 n, 2 n, 2 n} 

Alg 

{2 n + 1, 2 n + 1, 2 n + 1} 

a 2u 

{2n + 1, 2 n + 1, 2m }p 

T2 9 

(2 n, 2 n, 2m + l}p 

Tiu 

(2 n, 2 n, 2m}p 

Alg, Eg 

(2 n + 1, 2 n + 1, 2m + l}p 

A 2u i E u 

(2 n, 2m, 2k}p 

Alg, A 2 g, 

{2n + 1, 2m + 1, 2k + l}p 

Aiu, A 2u , 

(2n, 2m, 2k + l}p 

Tiu, T 2u 

(2n + 1, 2m + 1, 2k}p 

Tig, T 2g 


The anharmonic part of the potential reduces the symmetry of the system 
and the degeneracy of the energy levels is consequently diminished when 
A > 0 causing a splitting of the energy levels. For the lowest ones this 
splitting is given by 


V - 0 ^ A.lg 

v — 1 —> Tiu 

V — 2 —> Alg, Eg, T 2g 

v — 3 —y A 2u , Ti u , Ti u , T 2u 

V — 4 — y A\g, A\g, Eg, Eg, T\g, T 2 g, T 2 g. (5) 


This equation tells us that the lowest energy level is nondegenerate, the first 
excited energy level remains three-fold degenerate, the second excited energy 
level splits into a singlet a doublet and a triplet, the third excited energy level 
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splits into a singlet and three triples (two of them of the same symmetry Tin), 
etc. Note the alternating parity (either g or u) given by (—l) 1 ". 

For simplicity, in this paper we restrict ourselves to the case (3 = 0. This 


choice reduces the number of parameters in the potential-energy function but 


alters neither the symmetry of the problem nor the main conclusions drawn 
from it. 

3. Perturbation theory 

The purpose of this section is merely to carry out a simple calculation 
based on perturbation theory in order to verify the splitting of the energy 
levels outlined by equation (J5]). There are several strategies for obtaining 
the perturbation corrections when H 0 exhibits degenerate states. Here we 
simply obtain the first two perturbation corrections by inserting the trun¬ 
cated expansion E = + E^ A + E^ A 2 into the characteristic polynomial 

given by the secular determinant \H — EI\ = 0 and then solving the resulting 
equation for E ^ and E^ 2 \ This apparently impractical brute-force approach 
is sufficient for present purposes; the only subtlety being that each irrep is 
treated separately as discussed in section [4] (in particular with the basis set 
of subsection 14.3p . For the first energy levels we obtain the following results 



( 6 ) 


( 7 ) 


E\E g - 7 + -A — -A 2 + 0(A 3 ), 
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( 8 ) 




( 9 ) 


which are consistent with the prediction of group theory shown in equation 
(151) . It is worth noting that the degeneracy of the pairs of energy levels 
(E 2 A lg , and (Eia 2u , E 3 t 1u ) breaks at second order (the correction of 

first order being identical). However, such states offer no difficulty because 
we can treat each member of the pair separate from the other because they 
belong to different symmetry species. This strategy is one of the advantages 
of a possible systematic application of group theory to perturbation theory. 

4. Rayleigh-Ritz variational method 

This approach is based on a variational ansatz given by a finite linear 
combination of functions of a basis set B = {fo, fi, ■■■}'■ 


N -1 


i’ = Yi c ‘f‘- 


( 10 ) 


The variational coefficients Cj are chosen so that the variational integral 

F M H w 

is a minimum. The condition dE/dcj = 0 leads to the secular equations 


(fj \ H — E\ip) = 0, j 



( 12 ) 
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from which we obtain the coefficients c, : and approximate energies Ej(N), j = 
0,1,..., — 1. This result is valid for orthogonal as well as non-orthogonal 

basis functions fj, provided that in the latter case the functions are linearly 
independent. This approach always yields increasingly tighter upper bounds 
because Ej(N ) > Ej(N + 1 ) 


If the system exhibits symmetry it is convenient to choose a suitable 
basis set B s = {/jf, /f,...} for each irrep S. In order to obtain it we apply 
the projection operator P s (see the Appendix A for more details) to every 
fi E B and then remove the linearly dependent functions from the resulting 
set {P s fi, 2 = 0,1,...]-. One advantage of using the basis sets B s is that 
the dimension of the secular equations for a given accuracy is considerably 
smaller. Thus, instead of equations (TTUT) and (1T2T) we have 


N S -1 

v> s = E 

i =0 


and 


(f?\H-E s \iP s ) = 0, j — 0,1, , N s — 1 , 


(13) 


(14) 


for each irrep S. In this case we obtain upper bounds for each irrep exactly 
in the same way as before: Ej(Ns) > Ej(Ns + 1). 


4-1. Basis set of eigenfunctions of H 0 

One of the most convenient basis sets consists of linear combinations 
of eigenfunctions of H 0 adapted to the symmetry of H. In this case the 
variational method reduces to diagonalizing the matrix representation of the 
Hamiltonian H 5 in each orthonormal basis set B s = {<^f mn } adapted to the 
corresponding symmetry species S. We thus obtain approximate eigenvalues 




E s and eigenfunctions ip s that are linear combinations of the form 

C kmn<Pkmn ■ ( 15 ) 

k,m,n 

The coefficients cf mn and the approximate eigenvalues E s are given by the 
eigenvectors and eigenvalues of H 5 , respectively, and are also solutions of 
a secular equation similar to (fl4|) . The projection operators enable us to 
construct the symmetry-adapted basis sets B s in the following way 

P S Vkmn = U kmnVk'm'n'- ( 16 ) 

k' ,m' ,n' 

Since [H 0 , P s ] = 0 then v = 6 unless k + m + n = k! + w! + n' — v. 

After removing the linearly dependent functions and orthonormalizing the 
remaining ones the basis sets adapted to the symmetry of the Hamiltonian 
operator ([Tj) result to be 


A 


1 9 


iSP2n 2m 2m + ( f2m2n2m + i f2m2m2n) 


V^2n 2 n 2 n 

1 

7 ! 

i , 

\^2k2m2n + i P2n2k2m + ^>2m2n2k + ^>2m2k2n + i P2n2m2k 
+ ( / 9 2fc2n2?n) , 


(17) 


A 


2 g 


{ l f2k2m2n + ( f2n2k2m + <fi2m2n2k ~ <fi2m.2k2n ~ i f2n2m2k 


Vq 

^P2k 2 n 2m) j 


( 18 ) 


E n 


9 


A 


1 u 


1 1 

2m 2m P2m.2n2m l p2m2m.2n) i (^2m2n2m l p2m2m2n) 

(2</?2fc2m2n P2n2k2m P2m2n2k ) ) (V^2n2fc2m P2m2n2k) 

(2<y?2»n 2k 2n P2n2m2k P2k2n2m ) j ( fp2n2m2k P2k2n2m ) J(^9) 


T, 


19 


(9 7 2fc2m+12n+l — 9 7 2fc 2n+l 2m+l) , “^ (</?2m+l 2jfc 2ra+l — </?2n+l 2jfc 2ro+l) , 

( 20 ) 


( < / ? 2m+12n+12jt — 9 7 2n+l 2m+l 2fc) f * 




2 9 


\P2k2m+12m+l, P2m+12k2m+l, ( P2m+12m+12k} 

(P2k2m+12n+l + P2k2n+12m+l) i (P2m+12k2n+l + P2n+12k2m+l) , 


V2 


( ( P2m+12n+12k + P>2n+l 2m+l 2k) \ , 


( 21 ) 


—j= (P2k+12m+12n+l + ( P2n+12k+12m+l + P>2m+l2n+12k+l ~ P2m+12k+12n+l 

V6 


— P2n+12m+12k+l ~ P2k+12n+12m+l) , 


( 22 ) 


^2« • 

P2n+12n+12n+l 

—j= (ty?2fc+12m+12n+l + l p2n+12k+12m+\ + 9 7 2m+l 2n+l 2fc+l + <^2m+l 2fe+l 2n+l 

V 6 

+9 7 2n+12m+12fc+l + ^fc+l 2n+l 2m+l) j (23) 
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E u ■ 


y/6 


(2<^2n+l 2m+l 2m+l ^2m+l 2n+l 2m+l ^2m+l 2m+l 2n+l) > 


v/2 


(^2m+l 2n+l 2m+l 9^2m+l 2m+l 2n+l) 


V6 


(2y?2fc+12m+12n+l ~ ty?2ra+l 2fc+l 2m+l — ( f2m+l 2n+\ 2fc+l) j 


V2 


(9 7 2n+12fc+12m+l — ( / ? 2m+l 2n+l 2fc+l) f , 


(2</?2m+12fc+12n+l — </?2ra+l 2m+l 2fc+l ~ i f2k+l2n+l2m+l) , 


{tP2n+12m+12k+l ~ ty^fc+l 2n+l 2m+l) \ , 


(24) 


T, 




l9^2fc+12m2m> ^2m2fe+12m) 9^2m2m2fc+lj 

(ty?2fc+12m2n + l f2k+12n2m) , ( ^2m2k+l2n + f2n2k+12m) , 


—7= {^2m2n2k+\ + ty?2n2m2fc+l) 

V 2 


(25) 


T; 


2w 


(ty?2fc+12m2n — ^2fe+12n2m) , (</?2m 2fc+l 2ra — ^2n2k+12m) , 

~^( ( P2m2n2k+l ~ ( P2n2m2k+l)\ ■ (26) 

By convention tpijk means that all the subscripts are different, equal sub¬ 
scripts are indicated explicitly as, for example, in ipijj. 


4-2. Krylov space 

A particular basis set that spans what is commonly called the Krylov 
space is given by f t = H l f, where / is a suitably chosen function. If follows 


11 



from the properties of the projection operators discussed in the |Appendix A 
that P s ft = H l P s f = H l f s so that by simply choosing a seed function 
f s with the correct symmetry then the resulting basis set is automatically 
adapted to the corresponding irrep. We thus have B s = {ff — H l f s , i = 0,1,. 
for each irrep S and solve secular equations similar to (11411 . 

Suitable seed functions are 




f A 19 = exp [-a (x 2 + y 2 + z 2 )] , 


(27) 


f A29 = ( x 4 y 2 — x 4 z 2 — x 2 y i + x 2 z A + y^z 2 — y 2 z 4 ) exp [—a ( x 2 + y 2 + z 2 )] , 

(28) 

( 2z 2 — x 2 — y 2 ) exp [—a (x 2 + y 2 + z 2 )\ 

(x 2 — y 2 ) exp [—a ( x 2 + y 2 + z 2 )\ 


f Ea = 


(29) 


f Tl9 = 


(. xy 3 — x 3 y) exp [—a ( x 2 + y 2 + z 2 )\ 
(. xz 3 — x 3 z) exp [—a ( x 2 + y 2 + z 2 )\ 
(; yz 3 — y 3 z) exp [—a ( x 2 + y 2 + z 2 )\ 


(30) 


f T2g = { 


(31) 


xy exp [—a ( x 2 + y 2 + z 2 )\ 
xz exp [—a (x 2 + y 2 + z 2 )] , 
yzex p [—a ( x 2 + y 2 + z 2 )\ 
f Alu = ( x 5 yz 3 - x 5 y 3 z + x 3 y 5 z - x 3 yz 5 - xy 5 z 3 + xy 3 z 5 ) exp [-a (x 2 + y 2 + z 2 )\ , 

(32) 

fMu _ X y Zex p |^_ a ( x 2 + y 2 + z 2 )~\ , (33) 

xyz (2 z 2 — x 2 — y 2 ) exp [—a (x 2 + y 2 + z 2 )\ 
xyz (x 2 — y 2 ) exp [—a (x 2 + y 2 + x 2 )] 


/^ = 


(34) 


12 







( 35 ) 




xexp 

[—a ( 

x 2 + y 2 

+ ^~ 2 )] 


f Tlu = < 

yexp 

-a( 

x 2 + y 2 

+ ^)] 




zexp 

K 

-o( 

x 2 + y 2 

+ ^ 2 )] 



x ( y 2 — z 2 \ 

exp 

[—a ( x 2 

+ y 2 + 

* 2 )] 

= < 

y( 

z 2 - x 2 ) 

exp 

[—a ( x 2 

+ y 2 + 

* 2 )] 


> z ( 

x 2 — y 2 ] 

exp 

[—a ( x 2 

+ y 2 + 

^ 2 )] 


(36) 


where a > 0 is a nonlinear variational parameter that enables us to improve 
the accuracy of the results. 

We carried out a set of calculations with a = 1 and encountered a 
most surprising difficulty. For some reason (unknown to us at present) the 
Rayleigh-Ritz variational method in the Krylov space yields only one of the 
two T lu state triplets stemming from E (A = 0) = 7 and only one of the T 2g 
state triplets stemming from E (A = 0) = 11. We do not investigate this 
fact any further in this paper and just mention it in passing. This fact is 
surprising because we applied this approach in the past to other anharmonic 
oscillators and faced no such problem 7l|, 1721 ]. 


4■ 3. Non-orthogonal basis set 


In addition to the orthonormal basis set (jl 71126(1 we can also try a closely 
related symmetry-adapted non-orthogonal basis set given by 

B s = {P s x k y m z n exp [-a ( x 2 + y 2 + z 2 )] , k,m,n = 0,1,...} , (37) 

where a is a variational parameter that we set equal to unity for simplicity. 
In order to obtain a suitable basis set B s it is necessary to remove all the 
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linearly dependent functions produced by the application of the projection 
operator. Since this approach is quite straightforward for programming we 
chose it for present calculations of basis set dimensions N Al!l = 41, Ne 9 = 54, 
N Tlg = 66, N T2g = 102, N Mu = 16, N Tlu = 23, N T2u = 39. The first 
energy levels are shown in Figure |T] for 0 < A < 1. Present numerical 
results are accurate enough for the purpose of illustrating the splitting of the 
degenerate energy levels of the harmonic oscillator as A increases. They are 
consistent with both the splitting predicted by PGS ([5]) and the analytical 
perturbation results (|M9l). In fact, we obtained those perturbation results 
from the secular determinants obtained with the symmetry-adapted basis 
sets (137|) with a = 1/2 that is the exact value of this parameter for A = 0. 
This strategy was already described in Section [3J 


5. Conclusions 


Throughout this paper we have tried to illustrate the application of group 
theory to a quartic anharmonic oscillator with the symmetry Oh of the cube. 
We chose this particular example because of its great symmetry and also 
because it was treated before by means of a simpler symmetry-based ap¬ 


proach consisting only of parity and coordinate-permutation operators 


59]. 


The present application of group theory resorts to the 48 symmetry opera¬ 
tions shown in the table of characters in Table [Q and enables a systematic 
classification of the states of the oscillator in terms of the corresponding 
irreps. 

Group theory enables us to predict the splitting of the energy levels of the 
harmonic oscillator as the perturbation parameter A increases. This predic- 
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Table 1: Character table for group Oh 


0 h 

E 

8C 3 

6C 2 

6 C 4 

3C 2 (= Cl) 

i 

6S 4 

8S 6 

3o-fc 

6cr d 



A lg 

1 

1 

1 

1 

1 

i 

1 

1 

1 

1 


x 2 + y 2 + z 2 

to 

1 

1 

-1 

-1 

1 

i 

-1 

1 

1 

-1 



E g 

2 

-1 

0 

0 

2 

2 

0 

-1 

2 

0 


(2 z 2 — x 2 — y 2 , x 2 — y 2 ) 

Ti 9 

3 

0 

-1 

1 

-1 

3 

1 

0 

-1 

-1 

(Rx, Ry, Rz) 


T 2g 

3 

0 

1 

-1 

-1 

3 

-1 

0 

-1 

1 

( xz , yz , xy ) 


Aiu 

1 

1 

1 

1 

1 

-1 

-1 

-1 

-1 

-1 



Aiu 

1 

1 

-1 

-1 

1 

-1 

1 

-1 

-1 

1 



E u 

2 

-1 

0 

0 

2 

-2 

0 

1 

-2 

0 



Tin 

3 

0 

-1 

1 

-1 

-3 

-1 

0 

1 

1 

( x , y , z ) 


Tlu 

3 

0 

1 

-1 

-1 

-3 

1 

0 

1 

-1 




tion is verified by the actual calculation of the perturbation corrections. An 
advantage of the application of group theory is that we treat each symmetry 
subspace independently of the others. In this way both the calculation effort 
as well as the undesirably effect of the degeneracy are considerably reduced. 
Note, for example that the degenerate energy levels with v — 1 and v = 2 
can be treated as nondegenerate when one considers every symmetry species 
separately. 

Group theory also enables us to reduce the dimension of the secular equa¬ 
tions in the application of the Rayleigh-Ritz variational method. In addition 
to it, one obtains the eigenfunctions that are bases for the irreps. This 
fact facilitates, for example, the calculation of matrix elements of the form 
Onm = (V’nl ^ | V’m) f° r an y observable O. Given the symmetry of O one 


knows beforehand whether the matrix element O’P'T, is 


zero 


801421 . 


ysis leads, for example, to the selection rules for molecular spectra 


'his anal- 


8fM82j. 
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Figure 1: First eigenvalues of the Hamiltonian operator H = Px + P% + pi + x 2 + y 2 + 
z 2 + A ( x 2 y 2 + x 2 z 2 + y 2 z 2 ) with symmetry A\ g (solid line), T\ u (filled circles), E g (empty 
circles) , T 2g (filled squares), A 2u (dashed line), T 2u (empty squares), T lg (crosses). 


Appendix A. Simplified outline of Group Theory 

In what follows we develop a simplified and abbreviated version of some 
of the elements of group theory that are necessary for the present paper. 
A rigorous account of group theory is available in any of the books on the 


subject 8CH82] 


Here we are interested in the analysis of the symmetry properties of a 
physical system that are related to all the unitary operators Xj = X~ L that 
leave the Hamiltonian operator H invariant 

X t HX\ = H, (A.l) 

and assume that this set is finite 


G = {X u X 2t ...,X h } 


(A.2) 
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The product (composition) of these unitary operators is associative: (X t Xj)X k 
Xi(XjXk). The identity operator Xi, which satisfies 

Xrfi = X,X, = Xi, (A.3) 


leaves the Hamiltonian invariant and, consequently, belongs to G. It follows 
from (1A.1D that if Xi belongs to G then belongs to G too. The product 
of two operators also belongs to G as follows from 


X i XjH(X i X j ) Jf = XiXjHX^Xj = X, IIX- = H. 


(A.4) 


Because of all these mathematical properties the set G is a finite group. 

Two such operators ( or group elements) Xj and X k are said to be con¬ 
jugate if 

XiXjXt = X k , (A.5) 

for some A,- e G. If Xi and X k are conjugate to Xj then they are conjugate 
to each other. All the mutually conjugated elements of a group are collected 
into a class. 

We can construct a matrix representation X, of every operator Xi in 
terms of a basis 

B = (A.6) 

in the usual way 

Xifj = Y ( X 0/y fk, i = l,2,...,h. (A.7) 

k 

Of particular interest are the irreducible representations (irreps) obtained in 
terms of suitable basis sets 


B a 



a = 1, 2,..., m, 


(A.8) 
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such that 


*i/“=yh x .%/, 


kj J k "> 


j = 1,2, ...,l 0 


1,2,... ,h, a = 1,2,... ,m. 


k =1 

(A.9) 

The characters of the irreps are the traces of the corresponding matrix rep¬ 
resentations 


X? = tr (X“ 


(A.10) 


It can be proved that the number m of irreps equals the number of classes 
of group elements. 

In order to obtain a basis function f a for a given irrep we apply the 
corresponding projection operator 

P ° = \ £ (A?)* Xi, (A.11) 

i = 1 

to an arbitrary function f: f a = P a f. A projection operator P satisfies the 
following properties: P' = P and P 2 = P ; more precisely, any operator that 
satisfies these two properties is a projection operator. It follows from these 
properties that (P a f\ P a f) = (f \ P a f) < (/| /). 

In order to apply the equations above we should know how to express the 
effect of a group operator X t on a function /(x) of the cartesian coordinates 
x = ( x,y,z ). Rotations, reflections, etc can be expressed in matrix form as 
x' = Mx so that there is a one to one correspondence between the operators 
X t and the corresponding transformation matrices M,. If we write 


*i/(x) = / (M“ 1 x) , (A.12) 

then we have the mappings 


{Xi,Xj) (Mj, Mj) 
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X l X ] ->■ M,M j. 


(A.13) 


It follows from (jA.lj) that [H, X,] = 0 and, according to (1 A. 11 [1 . we con¬ 
clude that [H, P a ] = 0. 

In some cases it may be necessary to carry out an equivalent transforma¬ 
tion of the Hamiltonian of the system 


UHU ] = H, [A = U~\ 

The point group G for H 

G = {x l5 X 2 , ...,X h }, X t = UXirf, 
is isomorphic with G as follows from 

XiXj = UXiU^UXjU* = UXiXjW = XXj. 


(A.14) 


(A.15) 


(A.16) 


Table [T] shows the character table of the point group Oh- The first row 
shows the symmetry operations grouped into classes. E is the identity oper¬ 
ation, C n denotes a rotation by an angle 2tt/ n, i is the inversion operation, 
S n denotes a rotation by an angle 2tt jn followed by a reflection with respect 
to a plane perpendicular to the rotation axis and a indicates a reflection 
plane. The first column displays the irreps; those labelled by A, E, and T 
are one-, two- and three-dimensional (that is to say, l a = 1,2,3), respec¬ 
tively. The numbers are the characters xt an< 4 the last two columns show 
some basis functions for the irreps. This table summarizes some of the rel¬ 
evant ingredients for the construction of the projection operators P a shown 
in equation (lA.lljh The matrix representation of the elements of this group, 
which is necessary for the construction of the symmetry operations (1A.12D 
and projection operators flA.llH . is available elsewhere 73]. 
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